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ABSTRACT 

The single-particle Green's function of an interacting Fermi system with dom- 
inant forward scattering is calculated by decoupling the interaction by means of 
a Hubbard-Stratonowich transformation involving a bosonic auxiliary field <f> a . 
We obtain a higher dimensional generalization of the well-known one-dimensional 
bosonization result for the Green's function by first calculating the Green's func- 
tion for a fixed configuration of the (/> a -field and then averaging the resulting ex- 
pression with respect to the probability distribution V{(f) a } oc exp[—S e ff{<f> a }], 
where S e ff is the effective action of the <f> a -fve\A. We emphasize the approx- 
imations inherent in the higher-dimensional bosonization approach and clarify its 
relation with diagrammatic perturbation theory. 

1. Introduction 

The calculation of the single-particle Green's function G(k, u) of an interacting 
Fermi system in dimensions d > 1 is a very difficult problem, which can only be solved 
approximately. In conventional perturbative many-body theory G(k, u) is calculated 
by hrst expanding the irreducible self-energy E(k,cu) to some order in the effective 
interaction, and then re-summing the perturbation series by solving the Dyson equa- 
tion 1 . In some cases, however, the perturbation series for E(k,cu) is plagued by diver- 
gencies, which can only be cured if infinite orders in the interaction are re-summed. 
A well known example are interacting electrons in one spatial dimension. Fortunately, 
in d = 1 there exist several non-perturbative methods. Besides the Bethe-Ansatz 2 
and renormalization group methods 3 , the bosonization approach has been used with 
great success in one dimension 4 . Anderson's suggestion 5 that the normal-state proper- 
ties of the high-temperature superconductors are a manifestation of non- Fermi liquid 
behavior in d > 1 have revived the interest to develop non-perturbative methods in 
higher dimensions. To the best of our knowledge, the hrst attempt to generalize the 
bosonization approach to arbitrary d was due to Luther 6 . However, Luther's ideas did 
not receive much attention until Haldane 7 generalized them in such a way that boson- 
ization in d > 1 became a practically useful method for calculating the single-particle 
Green's function of a large class of interacting Fermi systems. Houghton et al. 8 as well 
as Castro-Neto and Fradkin 9 developed the bosonization approach further and applied 
it to some problems of physical interest. More recently, we have given an alternative 
formulation of higher dimensional bosonization 10 ' 11 , which is based on a generalization 
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of the functional bosonization approach developed in d = 1 by Fogedby 12 , and later 
by Lee and Chen 13 . A similar functional integral approach, which emphasizes more 
the mathematical aspects of bosonization in d > 1, has been developed independently 
by Frohlich et al. 14 . 

The advantage of the functional integral approach is that it can be used as a basis 
to calculate systematic corrections to the non-interacting boson approximation, which 
exist in any dimension if the energy dispersion is not linearized. Furthermore, in 
d > 1 so called "round-the-corner" processes (this terminology will become evident 
below) give rise to additional corrections to the leading bosonization result, so that 
in d > 1 bosonization is certainly not exact. In this paper we shall describe in detail 
the derivation of the single-particle Green's function by means of our functional bo- 
sonization approach. The approximations inherent in higher dimensional bosonization 
will be clarified and their validity will be critically discussed. We shall also derive 
the precise connection between the bosonization result for the Green's function and 
diagrammatic perturbation theory for the irreducible self-energy. In this way we shall 
recover a Ward-identity which has first been derived by Castellani, Di Castro, and 
Metzner 15 by means of a very different method. 



2. The Hubbard-Stratonowich transformation 



The crucial first step in the bosonization approach is the subdivision of the degrees 
of freedom close to the Fermi surface into sufficiently small boxes K% A of radial high A 
and cross section A d_1 . The index a labels the boxes in some convenient ordering. The 
cutoff A should be chosen sufficiently small so that within a given box the curvature of 
the Fermi surface can be locally neglected 7-11 . Given this partitioning of the degrees 
of freedom close to the Fermi surface, a general two-body interaction can be described 
by the Hamiltonian 

H mt = — J2 J2 c' : p\K ■ > C 1 ) 

Q aa' 

where : . . . : denotes normal ordering, and the /™ are coarse grained Landau interac- 
tion parameters. Here = J2k ©"(^OV'kV'k+q are the Fourier components of the local 
density operators associated with the boxes, where ip^ is the annihilation operator of 
an electron with momentum k. The cutoff function a (k) is unity for wave- vectors 
k G Kjl \ and vanishes otherwise. For simplicity, we shall consider spinless fermions. 
In an Euclidean functional integral approach, the exact Green's function of the system 
can then be written as 16 



G(k) = G(k,iCu n ) = -fj 



JD{0}e- s ™»'W ' v ' 

where the Euclidean action S mat {tp} is the following functional of the Grassmann 
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field ip, 



S ma tW = S {lp} + S lnt {tp} , 

s {4>} = /?E [~ l ^n + 6c] i>\i>k 



2V 



\ ^ \ ^ raa' a a' 



(3) 
(4) 

(5) 



q aa' 



Here <5c = Qc — is the energy dispersion measured relative to the chemical potential 
/i, and the local density operator is now represented by a composite Grassmann field 
Pq = J2k (k)4>k' l l>k+q- Throughout this work we shall use the convention that k = 
[k, iuj n ] and q = [q, iu m ], where the fermionic frequencies are Co n = 2ir(n + and 
the bosonic ones are u m = 2irm/f3. Here (3 is the inverse temperature, and V is the 
volume of the system. 

Defining the dimensionless parameters / q ""' = |/™', the interaction can be de- 
coupled by means of the following Hubbard-Stratonowich transformation over a bo- 
sonic auxiliary field c/>^, 



-S l nt{4>} 



exp 



<? aa' 



P- q P a q 



J X> {</>«} exp 


-\ E g 'rV-^' " * E g E„ 4>- q P q 


JV{<t> a } exp 







(6) 



Here / is a matrix in the patch indices, with matrix elements given by [/ ] aQ: ' = 

/™ . Eq.(6) is easily proved by shifting the c/^-field in the enumerator according to 

<f> q ~> 4> q -i £«'[/,]"" ' P"q and usin § the fact that [/,]""' = which is a trivial 

consequence of the hermiticity of the Hamiltonian. With the help of Eq.(6) we obtain 
from Eq.(2) 



G(k) = -(3 



/D{V>}D{^} 



where the decoupled action is given by 

sty, r} = So w + s 1 w>, r} + s 2 {<n 



(7) 



f8) 



wi 



ith 



q a 



1 



a 

q T-q 



W*} = 2 E EL£, 

<? aa' 



aa ja ja 



(9) 

(10) 
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Thus, the fermionic two-body interaction has disappeared. Instead, we have the prob- 
lem of a dynamic bosonic field (jf that is linearly coupled to the fermionic density. 
This field mediates the interaction between the fermions in the sense that integration 
over the c/^-field (i.e. "undoing" the Hubbard-Stratonowich transformation) generates 
an effective fermionic two-body interaction. In fact, because all interactions in nature 
can be viewed as the result of the exchange of some sort of particles, it is more general 
and fundamental to define the problem of interacting fermions in this way. This point 
of view has already been emphasized by Feynman and Hibbs 17 . 

The functional integral approach gives us the freedom of performing the fermionic 
integration before integrating over the c/^-field. To eliminate the fermions, we write 

SoW + s 1 {^,r} = -(3J2^l[G- 1 } kk ^ , (11) 

kk' 

where G~ l is an infinite matrix in momentum and frequency space, with matrix ele- 
ments given by the formal Dyson equation 

[G l ]kk' = [Go 1 ]^' ~~ \Y]kk' , (12) 
where Go is the non-interacting Matsubara Green's function matrix, 

[G ]kk' = hk'G (k) , G (k) = — — , (13) 

tu n - £ k 

and the generalized self-energy matrix is 

[vi**' = E© a ( k )H-*' , y q a = l s . (14) 

a P 

Recall that k denotes wave- vector and frequency, so that 8kk' = ^kk'^nn'- Choosing the 
normalization of the integration measure T>{tp} suitably, we have then 

J = detG -l = e TrlnG-i = g Xr In G^ 1 ^Tr \n[l — G V] ^ ^ 

and 

-fj [ V{^}^\e- S ^- S ^^ = [GUe^o 1 e Tr H i-GoV] (16) 



Hence, after integrating over the fermions the exact interacting Green's function in 
Eq.(7) can be written as an average of the diagonal element 

G(k) = f V{r}V{r}[G]kk = ([G] kk ) . (17) 



The normalized probability distribution V{(j) a } i 



is 



e Seff{^} 

" OT = /p»} t -w • (18) 
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where the effective action for the (/("-field contains in addition to the action £2 {</>"} 
given in Eq.(lO) a contribution due to the fermion determinant, 

Seff{<t> a } = S 2 {<f> a } + S kin {<t> a } , (19) 

with 

S km {r} = -Trln[l-G V] (20) 

Note that in Eq.(17) one first calculates the Green's function for a frozen configuration 
of the (/("-field, and then averages the resulting expression over all configurations of this 
field, with probability distribution given in Eq.(18). The above transformations are 
exact. Of course, in praxis it is impossible to calculate the exact interacting Green's 
function from Eq.(17), because (a) the matrix G~ l cannot be inverted exactly, (b) the 
kinetic energy contribution Ski n {<f> a } to the effective action of the (/("-field can only be 
calculated perturbatively, and (c) the probability distribution V{(j) a } in Eq.(18) is not 
Gaussian, so that the averaging procedure cannot be carried out exactly. The amazing 
fact is now that there exists a physically interesting limit where the difficulties (a), (b) 
and (c) can all be overcome. The above method leads then to a new non-perturbative 
approach to the fermionic many-body problem. The highly non-perturbative character 
of this approach is evident from the fact that in d = 1 the exact solution for the Green's 
function of the Tomonaga-Luttinger model 4 can be obtained in this way 13 ' 18 . In d > 1, 
this methods leads to a straight-forward generalization of the bosonization approach to 
arbitrary dimensions. In the next section we shall discuss in detail how this calculation 
is carried out in praxis. 

3. Calculation of the Green's function 

3.1. The Gaussian probability distribution 

In general the above kinetic-energy contribution to the effective action can only be 
calculated perturbatively by expanding 

00-1 00 
S km {cf> a } = ~Tr ln[l - G V] = £ -Tr [6 v] n = £ S kin , n {<j> a } , (21) 

n=l n=l 

and truncating the expansion at some finite order. Within Gaussian approximation 
all terms with n > 3 in Eq.(21) are neglected, so that one approximates 

S k in{<t> a }*Tr[G V] + l -Tr[6 V] 2 . (22) 
Carrying out the traces, it is easy to show that within Gaussian approximation 

%{f}^EC^ + % 2 {f} , (23) 
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with 

SejjAr} = ^ e EtLCr' + K a 'mr- q ^ , (24) 

Z P q aa> 



where / is again a matrix in the patch-indices, with matrix elements given by the 
Landau parameters defined in Eq.(l). Here = J2k © a (k)/(^k) is the number 
of occupied states in box K% x , where f(E) = p je +1 is the Fermi function. The 
polarization part in Eq.(24) is given by 

ng a, (g) = -^E e"(k)0"'(k + q) { (a+q) ~ /(6c) . (25) 

For |q| <^ kp we may approximate Q a (k.)Q a> (k + q) ~ 8 aa 'Q a (k) } so that to leading 
order in |q|/&p we have in any dimension 

Il a /(q) « 6 aa 'lIZ(q) , W {q) = u a -^—^— , (26) 

V" • q - lU m 



where 



1 BN a 1 
1 uiy o _ _ qc 

V d/i ~ V ^ 



5/(6 



% (27) 

is the local density of states associated with patch a, and v" is the local Fermi velo- 
city. Note that the approximation in Eq.(26) is valid for small |q|/&p but for arbitrary 
frequencies, and that the local density of states v a is in general a cutoff-dependent 
quantity. In the work by Houghton et al. 8 and Castro Neto and Fradkin 9 it is im- 
plicitly assumed that the Gaussian approximation is justified. However, in none of 
these works the corrections to the Gaussian approximation are calculated, so that the 
small parameter which actually controls the accuracy of the Gaussian approximation 
remains hidden. Recently we have calculated this parameter by means of our func- 
tional bosonization approach 11 . Note that the Gaussian propagator of the c/^-field is 
simply given by the RPA-interaction matrix, 



where the elements of the matrix Ho(?) are given in Eq.(25). 

In the exactly solvable one-dimensional Tomonaga-Luttinger model 4 the bosonized 
Hamiltonian is known to be quadratic, so that in this case we have exactly 

-Tr\n[l - G V] = Tr [d V] +^Tr [d V] 2 . (29) 

All higher order terms vanish identically due to a large scale cancellation between 
self-energy and vertex corrections, which has been discovered by Dzyaloshinskii and 
Larkin 19 . A few years later T. Bohr gave a more readable proof of this cancellation 20 , 
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and formulated it in form of a theorem which he called the closed loop theorem. As 
shown in Ref. 11 , the cancellations responsible for the validity of Eq.(29) in d = 1 exist 
also in higher dimensions provided the following two approximations are made 21 : 

(Al): High density-limit or small momentum-transfer limit. 

The interaction should be dominated by momentum transfers |q| ^ q c <C kp. 
Because the fields (\) a q mediate the interaction, this condition is equivalent with the 
statement that the short wave-length Fourier components c/>^ with |q| ^ q c can be neg- 
lected. In this case we may ignore processes that transfer momentum between different 
patches, provided the size of the patches is larger than q c . Note that these "around- 
the-corner" processes exist for any d > 1, because there are always some neighboring 
patches which can be connected by arbitrarily small momentum transfers. Formally, 
the existence of the small parameter justifies the neglect of around-the-corner pro- 
cesses. Evidently, the Gaussian approximation would be completely uncontrolled for 
the Hubbard model, where fluctuations on all length scales are important. 

(A2): Local linearization of the energy dispersion at the Fermi-surface. 
Suppose we linearize the energy dispersion within a given box K% A by approxim- 
ating 

q = e ka+q -^«v a -q , (30) 

where k" is a vector on the Fermi surface that points to the center of patch a. Local 
linearization amounts to replacing Eq.(4) by 

5oW«)9EE 0a ( k )H^ + va -( k - ka M^ • ( 31 ) 

k a 

Note that in this approximation the Fermi surface is approximated by a collection of 
flat d — 1 dimensional hyper-surfaces, i.e. planes in d = 3 and straight lines in d = 2. 

If the approximations (Al) and (A2) are valid, then the generalized closed loop 
theorem discussed in Ref. 11 implies that the probability distribution V{(j) a } can be 
approximated by a Gaussian, 

POT »W-}= j-p W .-«.,„<,., • P2) 

At zero temperature the first term in Eq.(23) involving the zero wave- vector and 
frequency component of the c/^-field can be ignored for the calculation of correlation 
functions at finite wave- vectors or frequencies. 

3.2. The Green's function for fixed configuration of the (f) a -field 



7 



In order to calculate the Green's function from Eq.(17), we need to invert G~ l . 
We proceed in two steps. We first show that the condition (Al) means that G~ l 
is approximately block-diagonal, with diagonal blocks (G a )~ l labelled by the patch 
indices. Therefore the problem of inverting G~ l can be reduced to the problem of 
inverting each diagonal block separately. 

Block diagonalization. 

The quadratic form defining the matrix elements [G~ l ]kk' in Eq.(ll) can be written 

as 

S {^} + SMA a } = -PT.^l\G~%M^k +q , (33) 

kq 

with 

[Gl^^eikiju^-^j-y;] , (34) 

a 

where ^ is defined in Eq.(30). The cutoff function a (k) groups the matrix elements 
of the infinite matrix G~ l into rows labelled by the patch index a. To see this more 
clearly, consider a spherical Fermi surface in d = 2. As shown in Fig. I, we partition 
the degrees of freedom in the vicinity of the Fermi surface into n boxes and label 
neighboring boxes in increasing order. Evidently the group of matrix elements 




Fig. 1. Subdivision of the degrees of freedom close to spherical Fermi surface in d = 2 into boxes. 



associated with a given index a in Eq.(34) correspond to the horizontal stripes in 
the schematic representation of the matrix G~ l shown in Fig. 2(a). The width of 
the diagonal band with non-zero matrix elements is determined by the range q c of 
the interaction in momentum space, because the vanishing of the interaction f^ a for 
|q| £ q c implies that the held V q a mediating this interaction must also vanish. Noting 
that by assumption (Al) q c <C kp } we may choose the patch cutoffs A and A such that 

g c <A,A<A; F . (35) 
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- k' » k' 

1 2 n 




(a) ( (b) 

Fig. 2. (a) Schematic representation of the matrix G _1 defined in Eq.(34) for d = 2. Only the wave- 
vector index is shown, i.e. each matrix-element is an infinite matrix in frequency space. Regions 
with non-zero matrix elements are shaded, (b) Diagonal blocks and around- the corner processes 
(represented by black triangles). 



As shown in Fig. 2(b), in this way the matrix G~ l is subdivided into block-matrices 
associated with the patches such that G~ l is approximately block-diagonal. The two 
black triangles in the upper right and lower left corner of the matrix in Fig. 2(b) 
represent scattering processes between patches 1 and n. Because these patches are 
adjacent, they can be connected by arbitrarily small momentum transfers. The crucial 
approximation is now to neglect all matrix elements describing momentum transfer 
between different boxes, i.e. the around-the-corner processes. These are located in the 
black triangles of Fig. 2(b). The justification for this step is that the relative number 
of matrix elements representing such processes is small as long as the condition (35) 
is satisfied. In d > I dimensions the relative number of around-the-corner matrix 
elements for any given is of the order of A ^-i A <C 1. Note that this approximation 
makes only sense if the patch cutoffs are kept finite and large compared with the range 
of the interaction in momentum space! Once we have disposed of the matrix elements 
in the black triangles of Fig. 2(b), the matrix G~ l is a direct sum of diagonal blocks 
(G") -1 , a = 1, . . . , n, so that 

[(TV = E 0a ( k ) 0a ( k O[(GT 1 W , (36) 

a. 

where the matrix (G a )~ l is the diagonal block of G~ l associated with patch a, 

[{&)-%» = 8 kk \iu n -(£_ va ]-V£_ k , . (37) 

Thus, the problem of inverting G~ l is reduced to the problem of inverting each diagonal 
block separately. The diagonal elements of G are then simply given by 

[G^* = £ 0a ( k )[<H* • (38) 

a. 
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Note that G a is still an infinite matrix in frequency space, so that the quantum dy- 
namics is fully taken into account. 

Although the relative number of matrix elements describing around-the-corner pro- 
cesses is small, we have to make one important caveat: Possible non-perturbative 
effects that depend on the global topology of the Fermi surface cannot be described 
within this approximation. For example, in d = 2 each patch has two neighbors, but 
the first and the last patch are adjacent, so that there exist also around-the-corner 
processes between patch 1 and patch n, which give rise to the off-diagonal triangles in 
the lower left and upper right corners of Fig. 2. More generally, in higher dimensions 
the off-diagonal around-the-corner blocks are distributed in a complicated way over 
the matrix G~ l . The effect of these sparsely distributed around-the-corner blocks is 
difficult to estimate, and we are implicitly assuming that they do not lead to qualit- 
atively new effects. We would like to emphasize that this is an assumption which is 
implicitly also made in the operator bosonization approach 8 ' 9 , and in the Ward-identity 
approach by Castellani, Di Casto and Metzner 15 . 

Inversion of the diagonal blocks. 

Up to this point we have not linearized the energy dispersion, so that the above 
block diagonalization is valid for arbitrary dispersion £q. The crucial advantage of the 
subdivision of G~ l into blocks is that within a given block it is allowed to linearize 
the energy dispersion, ^ Ri v™ ■ q. This is the approximation (A2). By linearizing 
the energy dispersion, we ignore locally the curvature of the Fermi surface. This is 
accurate if the patches are sufficiently small, so that the variation of the direction of 
the local normal vector can be ignored within a given patch. Note, however, that the 
size of the patches should be chosen large enough to satisfy Eq.(35). Linearization 
is always justified if the curvature of the Fermi surface is intrinsically small, so that 
even for large |q| the corrections to ^ v" • q can be neglected. In the latter 
case it is sufficient to cover the Fermi surface with a small number of patches. For 
example, for quasi-one-dimensional chain-like materials the Fermi surface consists of 
two disconnected almost flat pieces. As shown in Ref. 22 , in this case the magnitude 
of A can be chosen of the order of kp. 

Once the linearization has been made, it is possible to invert the diagonal block 
(G a )~ l exactly. Note that G a is still an infinite matrix in frequency space. Shifting 
the wave-vector labels according to k = k" + q and k' = k" + q', the diagonal block 
G a is determined by the equation 

£ [hAGZWr 1 - VjL?] [G a ]^ = <W • (39) 

q' 

where Gq(c[) = [iuj n — v" • q] -1 (For simplicity we have introduced the collective label 
q = [q, iu n ].) The important point is now that Eq.(39) is first order and can be solved 
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exactly by means of a trivial generalization of a method due to Schwinger . Defining 
a (r,r',T,r') = £ e i ^ r -^ e -^'- r '-^[G a ] qql (40) 

I qq' 

V a (r,r) = ^yta-r-^V" , (41) 

it is easy to see that Eq.(39) is equivalent with 

[-d T + iv a - V r -V a (r,r)]0 a (r,r',T,<r') = 8(r - r')8*(r - t') , (42) 

where 

8*(t-t') = I^ e --"(— ') . (43) 

P n 

Note that the Fourier transformation in Eq.(40) involves fermionic Matsubara frequen- 
cies, because £? a (r, r', r, r') should be anti-periodic in each imaginary time variable. In 
contrast, V a (r, r) should be a periodic function of r, so that V q a depends on bosonic 
Matsubara frequencies. We now substitute the ansatz 

g a {r, r', r, r') = G£(r - r', r - r ') e * a (^)-* a (rV) ( 44 ) 
into Eq.(42), where Gq (r — r', r — r') satisfies 

[-<9 T + zv a - V r ]G£(r-r',r-r') = (J(r - r')<T(r - r') . (45) 

This yields 

[-d T + jv" • V r - V a (r, r)] £T (r, r', r, r') = 

8(r - r')8*(r - t') + g a (r, r', r, r') {[-d T + zv a • V r ] $ a (r, r) - V a (r, r)} . 

(46) 

Comparing Eq.(46) with Eq.(42), we see that our ansatz is consistent provided $ a (r, r) 
satisfies 

[-d T + iv a - V r ]$ a (r,r) = V a (r,T) . (47) 

Eqs.(45) and (47) are first order linear differential equations, which can be easily 
solved via Fourier transformation 

1 t(q-r-u„T) 

G «<-» = ^S^— ~ q • (48) 

i(q-r-c<; m T) 

$ a (r,r) = Y- V" . (49) 
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Having determined Gq(y,t) and $ a (r,r), the diagonal blocks (G a ) 1 are inverted, so 
that G a is explicitly known as functional of the c/^-field. 

3.3. Gaussian averaging: calculation of the Debye-Waller 

Combining Eqs.(17), (40), (44) and (49), and using the fact that averaging restores 
translational invariance in space and time, we conclude that the interacting Matsubara 
Greens-function is given by 

G(k) = V0 a (k) / dr / /? C /re- 8 [( k - ka )™ T ]G'^(r,r)(e $a ( r ' T )- $a ( ' )) . (50) 

„ J J0 X ' S eff,2 

Using Eqs.(14) and (49) we may write 

$«(r,r)-$«(0,0) = ^^(r,r)^ , (51) 



with 



l _ e -i(q-r-w m r) 



(52) 



(3 l iu m - v" • q 

The problem of calculating the interacting Greens-function is now reduced to a multi- 
dimensional Gaussian integration, which simply yields the usual Debye-Waller factor, 



D $ a (r,r)-$ a (0,0) 



1 



exp 



exp 



*/S, 



eff,2 



J°(r,r)J q a (r,r) 



P 



9V U^ PA rJ- q (r^)J q a (r,T) 

q 



(53) 



where we have used the fact that the Gaussian propagator of the c/^-field is accord- 
ing to Eq.(28) proportional to the RPA-interaction. For consistency, in Eq.(53) the 
polarization contribution to [f^ PA ] aa should be approximated by its leading long- 
wavelength limit given in Eq.(26), because in deriving Eq.(53) we have neglected 
around-the-corner processes. Using 



J"(r,T)J q a (r,T) 



2 1 — cos(q • r — cu m r) 

(3 2 (lLU ri 



we conclude that 



D $ a (r,r)-$ a (0,0) 



V" • q) 



Q a (r,r) 



where the Debye-Waller factor Q a (r } r) is given by 



(r,r) = R a - S a {r,T) , R a = lim5' a (r,r) , 

r,r— >0 



(54) 
(55) 

(56) 
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with 

1 _ f* PA ' a cos(q • r - u m r) 



S a (^r) = —Y J J -^-. V " " ' • (57) 



(3V Y fa™ ~ v" • q) 2 



Here f p ,a = [/ ] aa is the diagonal element of the RPA-interaction matrix. An 
important special case is a patch-independent bare interaction, i.e. [/ ] aQ: = / q for 
all a and a'. Then it is easy to see that f PPA,a can be identified with with the usual 
RPA-interaction, 



h -h -i + / q n ( g ) ' ^ J " /q ' ( j 

where Ho(q) = J2 a (?) 1S the total non-interacting polarization. 



In summary, we obtain for the full Matsubara Green's function of the interacting 
many body system 

G{k) = ^2 ® a (k)G a (k - k a , iu n ) , (59) 

a 

where 

G a (q) = G a (q,iCo n ) = J dr dTe- t{(l - r -^G a (r,T) , (60) 

with 

G a (r,T) = G a Q {v,T)e Qa{r ^ . (61) 

Shifting in Eq.(59) k = k" + q and choosing |q| small compared with the cutoffs A 
and A that determine the size of the box K£ x , it is easy to see that only the term 
a' = a in the sum contributes, so that (after renaming again a' — > a) 

G(k a + q } iu n ) = G a (q } iu n ) , |q|«A,A . (62) 



4. The connection with diagrammatic perturbation theory 

In this section we shall derive from Eqs.(59)-(61) an expression for the irreducible 
self-energy of the many-body system and compare it with the skeleton diagram. In 
this way we also elucidate the precise relation between higher dimensional bosonization 
and the Ward-identity approach by Castellani, Di Castro and Metzner (CCM) 15 . 

4-1. The integral equation 

Let us apply the differential operator — d T + iv a ■ V r to the bosonization result for 
the patch Green's function G a (r, r) in Eq.(61). Using the fact that the application of 
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this operator to Gq(y,t) generates the usual ^-function (see Eqs.(45) and (43)), it is 
easy to show that 

[-d T + iv a -V r + X a (r- r', r - r')] G a (r -r',r- t') = 8(r - t')8*(t - t') , (63) 
with 

X a (r-r',T-T') = -[-d T + iv a - V r ]Q a (r-r',r-r') . (64) 
From the explicit expression for Q a (r } r) given in Eqs.(56) and (57) we find 

i fRPA,a 

X«(r,r) = w j:^-^x; , X ; = 1 -b—. . ,65) 

In Fourier space Eq.(63) becomes 

[iu n - v" • q]G a (q) + -±- J2 X?_ 5 ,G a {q') = 1 , (66) 



or equivalently 



rRPA,a 

[id; n -v--q]G a (q,id; n ) = 1-^-5:^ '^'f ■ ( 67 ) 

pV iu n _ n , - v" • (q - q') 

Recall that q = [q, iu m ] involves bosonic Matsubara frequencies, whereas the label 
q = [q, iuj n ] depends on fermionic ones. Because the difference between two fermi- 
onic Matsubara frequencies is a bosonic one, the kernel X~_~, in Eq.(66) depends on 
bosonic frequencies. In the zero-temperature limit Eq.(67) is equivalent with the in- 
tegral equation given in Eq.(13) of the work by Castellani, Di Castro, and Metzner 15 . 
Our bosonization approach reduces the solution of Eq.(67) to the standard problem of 
solving a linear partial differential equation (Eq.(46)) and calculating a Debye- Waller 
factor in a Gaussian integral. In obtaining this solution, the non-trivial atlas with local 
coordinate systems on the Fermi surface has played an important role. This patch- 
ing construction is essential to exhibit the large-scale cancellation between self-energy 
and vertex corrections, which according to the closed loop theorem are guaranteed to 
happen in arbitrary d if the conditions (Al) and (A2) listed above are satisfied. An 
important difference between our method and the Ward-identity approach of CCM 
is that these authors do not make use of the patching construction and the associ- 
ated non-trivial atlas on the Fermi surface. Instead, CCM work with a single rigid 
coordinate system, with origin at the center of the Fermi sphere. In this coordinate 
system the integrations seem to be technically more cumbersome. Note also that the 
final expression for the Green's function derived by CCM explicitly depends on some 
ultraviolet cutoff, so that it is not obvious how for sufficiently long-ranged interactions 
the result becomes cutoff-independent. On the other hand, with the patching construc- 
tion one can take advantage of the fact that for interactions / q that are dominated by 
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|q| <C kp only the local curvature of the Fermi surface within a given patch determines 
the accuracy of the approximations. Therefore Eqs.(59)-(61) are accurate in arbitrary 
dimension provided the approximations (Al) and (A2) discussed above are justified. 
In this case our approach leads to cutoff-independent results for physical correlation 
functions 10 ' 11 ' 18 . 

J f .2. The Ward identity 

In diagrammatic perturbation theory it is sometimes convenient 24 to define so 
called skeleton diagrams in order to exhibit the structure of the perturbation series 
more clearly. The skeleton diagram for the exact self-energy is shown in Fig. 3. In the 




Fig. 3. Skeleton diagram for the irreducible self-energy. The thick wavy line denotes the exact 
screened effective interaction /* , the shaded triangle is the exact three-legged vertex, and the solid 
line is the exact Green's function. 



Matsubara formalism, this diagram represents the following expression, 

s (*0 = -jy £ fUMk; k - k')G(k') . (68) 

f 

The exact effective interaction f* is related to the bare interaction via f* = -4^7, where 

J q J q t (q1 7 

e(q) is the exact dielectric function. By definition, the vertex function A(k; q) is the 
sum of all diagrams with three external ends corresponding to two solid lines and a 
single interaction line. Because G(k') on the right-hand side of Eq.(68) depends again 
on via the Dyson equation, Eq.(68) is a complicated integral equation, which 

can only be solved approximately. Moreover, the formal kernel f£_ k ,A(k; k — k') of 
this integral equation is again a functional of G(k) } so that it cannot be calculated 
exactly unless the entire perturbation series has been summed. 

For better comparison with the self-energy calculated within our bosonization ap- 
proach, let us shift again k = k" + q and k' = k" + q', so that wave- vectors are 
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measured with respect to the local coordinate system a. Defining 

G(k a + q,iu n ) = G a (q) , (69) 
E(k" + q,z^) = E"(g) , (70) 
A(k a + q,iu n ;q-q',iu n _ n ,) = A a (q,q-q') , (71) 

the skeleton equation (68) reads 

x a (q) = -^^ri-^ a m-?)G a (?) , (72) 

while the Dyson equation can be written as 

[G a m~ l = [G^q)}- 1 - £"(<?) • (73) 

Let us now determine the skeleton parameters that correspond to our bosonization 
result for the Green's function. Starting point is the integral equation (66). Noting 
that after linearization iu n — v" • q = [Gq (q)]~ l and dividing both sides of Eq.(66) by 
G a (q) } we obtain 

1 X a 

[G a (q)Y l = [GZffl- 1 + ^ E ^ff-iq ) • (74) 

Comparing this with Eq.(73), we conclude that in our bosonization approach the self- 
energy satisfies 

™=-w$<g%rw ■ <75 » 

From Eqs.(72) and (75) we obtain 

va rRPA,a 

Hence, the approximations inherent in our bosonization approach amount to replacing 
the exact effective interaction /* by the RPA-inter action, and approximating the vertex 
function by 

A a (q; q - q ) = , (77) 

[iun-n' ~ v a • {q-q')]G a {q) 

which is equivalent to 

[ l u m ,-v a - q '}A a (q;q') = [G a (q)}- 1 . (78) 

The important point is that the left hand side of Eq.(78) depends again on the exact 
Green's function. Such a relation between a vertex function and a Green's function is 
called a Ward identity. Using the definition (69) and the fact that after linearization 
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Gq(c[) = — Gq(— q), it is easy to see that Eq.(78) is equivalent with the Ward identity 
derived by CCM 15 . Thus, although within the bosonization approach the dielectric 
function is approximated by the RPA-expression, bosonization does not simply re- 
produce the usual RPA self-energy, because it sums in addition infinitely many other 
diagrams by means of a non-trivial Ward-identity for the vertex function! The ana- 
lytic expressions for these diagrams can be easily obtained by iterating the integral 
equation (67). Let us recapitulate how we have obtained Eq.(78): Starting point was 
the bosonization result for the Green's function of the interacting many body system 
in Eq.(61). By simple differentiation we have obtained the integral equation (67). Fi- 
nally, we have shown that this integral equation implies the Ward identity (78). The 
strategy of CCM was to perform these steps in precisely the opposite order. 

5. Conclusions 

In this work we have shown that by simple Gaussian integration one can obtain a 
non-perturbative expression for the single-particle Green's function which sums infinite 
orders in perturbation theory in a consistent way. Our approach is controlled in the 
high density limit for arbitrary dimensions if the interaction is dominated by small 
momentum transfers. The fact that in one dimension we obtain the exact solution 
of the Tomonaga-Luttinger model shows that our approach does not assume a priori 
that the system is a Fermi liquid. Thus, with the method described in this paper it is 
possible to study strongly correlated Fermi liquids with small quasi-particle residue, 
as well as possible non Fermi liquid states in higher dimensions. For example, we 
have applied our functional bosonization approach to quasi-one-dimensional metals 22 , 
to electrons coupled to transverse gauge fields 25 , to strongly coupled electron-phonon 
systems 18 , and to the problem of electrons moving in a stochastic medium 18 . Although 
the derivation of Eqs.(59)-(61) was rather straight-forward, the evaluation of the full 
momentum- and frequency dependent spectral function from these expressions is a 
very difficult mathematical problem, which so far has not been completely solved. A 
detailed discussion of Eqs.(59)-(61) and applications to problems of physical interest 
will be published in a forthcoming book 18 . 
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